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Abstract 

A two-parameter family of asymmetric exclusion processes for particles on a one- 
dimensional lattice is defined. The two parameters of the model control the driving 
force and an effect which we call pushing, due to the fact that particles can push each 
other in this model. We show that this model is exactly solvable via the coordinate 
Bethe Ansatz and show that its iV-particle 5-matrix is factorizable. We also study 
the interplay of the above effects in determining various steady state and dynamical 
characteristics of the system. 
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1 Introduction 

The asymmetric simple exclusion process (ASEP) is the simplest and most studied model 
of interacting particle systems in one dimension (see [1-5] and references therein). This 
model can be related via suitable mappings to many different physical models ranging from 
interface growth [6,7] to problems of traffic flow [8-11]. In this paper, however we would 
like to look at it only as a model for a collection of random walkers interacting with each 
other by simple exclusion. Each particle hops with rate R(L) to it's right (left) neighboring 
site if this site is empty, otherwise it stops. When the hoping rates R and L are equal, one 
deals with the symmetric exclusion process. This case, models diffusion of particles on a one 
dimensional line in the absence of a driving force. The other extreme case (i.e. L=0), is 
usually called the totally asymmetric exclusion process (TASEP). 

The asymmetric exclusion process has been extensively studied in the past few years as a 
prototype of a one dimensional system far from equilibrium for which some exact results 
can be obtained. Among the many aspects of this problem which have been studied, one 
can mention the mean field solution of the steady state and the phase structure [12,13], the 
exact steady state of the open [14], and closed chain [15], the effect of impurities [16-18], 
exact calculation of some dynamical properties [19], the effect of different kinds of updatings 
[20-21], and finally the exact calculation of conditional probabilities via the coordinate Bethe 
Ansatz [22]. 

In ref. [23] we considered the TASEP and added a new element into this process, namely the 
possibility that a particle pushes the particles in front of it with a rate depending on the 
number of these particles. Being interested in exactly solvable models, we found that if a 
particle can push a collection of n-adjacent particles in its immediate neighborhood with a 
rate given by 

1 + - + h (- 

where 0</i=l — A < 1, then the problem allows an exact solution via the coordinate Bethe 
Ansatz. The fact that r n decreases with n is physically natural, although its functional form 
may seem peculiar. This particular form is dictated by our demand of the exact solvability of 
the model. Varying the parameter /z from to 1, one can then smoothly interpolates between 
the TASEP and the so called drop-push model [24]. 

Remark: Here, by an exact solution we mean an exact determination of the time depen- 
dence of n-particle conditional probabilities on an infinite lattice. We also show that the 
steady state of the system on a ring is one in which all of the configurations have equal 
weights. The steady state of the system on an open chain is not known at present. Inciden- 
tally, the technique of matrix product ansatz does not work for this latter problem since the 
Hamiltonian of the process is highly nonlocal. 
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In the present work we consider the partially asymmetric case and study the combined 
effect of pushing and driving. Again we demand exact solvability and follow the strat- 
egy of our previous work [23]. The basic objects we are interested in are the probabilities 
P(x\, X2, ■ ■ • i xn] i 1 2/1 , y%i • ■ ■ , Un) for finding at time t, particle 1 at point x%, particle 2 at 
point x 2 ,- • •, given their initial positions at points y%,y2, 4 --, respectively. For notational 
convenience, in the following we suppress the initial coordinates whenever no confusion can 
arise. The physical region for the coordinates is x\ < X2 < ... < xn- If one wants to take into 
account a type of pushing effect, one should write a large number of equations depending 
on which and how many of the particles are adjacent to each other. The number of cases 
and hence the number of equations grows rapidly as the number of particles increases. The 
basic idea is to write one single equation supplemented by a particular boundary condition 
and to see what kind of pushing effect emerges. The master equation is 

d 

—P(xx,X2,-",^N',t) = R[P(xi - l,x 2 , ■ ■ ■ ,x N ;t) H h P{xi, x 2 , ■ ■ • , x N - 1; t)\ 

+L[P(x 1 + 1, x 2 , ■ ■ ■ , x N ; t) H h P(xi,x 2 , ■ ■ ■ , x N + 1; t)] 

-NP{ Xl ,x 2 ,---,x N ;t). (2) 

In the following we re-scale time so that R + L — 1. We assume that this equation holds in 
the whole physical region. When any of the coordinates of the left hand side are adjacent 
(say x i+ i — Xi + 1) , points on the boundary of the physical region appear on the right hand 
side of (2). We fix the value of these terms by the following boundary condition: 

P(x, x; t) = XP(x, x + l;t) + fiP(x - 1, x; t), (3) 

where for simplicity we have suppressed all the other coordinates. This idea is most apparent 
in sectors of low particle number. For example, in the two particle sector, combination of 
(0) and (|3]) yields for adjacent particles 

d 

—P(x,x + 1) = R[P(x - l,x + 1) +/j,P(x - l,x) - (1 + n)P(x,x + 1)] 

+L[P(x,x + 2) + XP(x + l,x + 2)-(l + X)P(x,x + l)], (4) 

which means that a particle hops freely to the right and left with rates R and L and pushes 
its neighboring particle to the right and left with rates Rfi and LA, respectively. We will 
show that in the general case, equations (@) and (|[) imply that the following processes : 

(5) 

011---1, — >11_^_1.0 (6) 

n+1 n+1 
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occur respectively with rates r n and l n given by 



f'n — R . / \\ n ' V/ 

1 + - H h - 



and 



Remark: Interchanging particles and holes, this model can be seen to be equivalent to a 
model in which each particle does not only hops to its immediate neighboring sites but to 
any other vacant site, with rates depending on the hopping distance; i.e., in the equivalent 
model the following processes occur with rates r n and /„, respectively 

00^_01 — >100^_0 (9) 

n+1 n+1 

100^_0 — >00_^_01 (10) 

n+1 n+1 

. The basic parameters of the model are R = 1 — L, and fx — X — 1. As we increase fx from 
to 1, the effect of pushing to the right increases and that to the left decreases. For single 
particles the hopping rates to the right and left are still R and L respectively. Thus the 
driving is controlled by R or L and the pushing by fx or A. 
A few special cases are worth mentioning: 

When fx = A, the pushing effect to the right and left are equal and the rates are 

r n = R^— l n = L — ^— . (11) 

n+1 n+1 

In this case, the asymmetry is controlled only by driving. 

When fx — 1, then we have maximum pushing to the right and no pushing to the left. The 

rates are 

r n — R-i ^0 = L, l n> Q = 

Single particles can hop to the left, but can not push other particles to the left. 
When L = 0, particles hop to the right only, and by varying fx form to 1 we interpolate 
between TASEP and the drop-push model [24]. Another way to understand the difference 
of the two sources of asymmetry in this problem is to note that the effect of driving can 
at least at long times be completely removed by going to an appropriate frame of reference 
while that of pushing can not. To see this, consider the transformation 

Xi -> x\ := Xi - Vt, (12) 

which is a Galilean boost. Actually, this is not an allowed transformation of our problem, 
since the x^'s are integers, whereas t is real. However, if the probability distribution is 
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sufficiently slowly varying (e.g. at long times), then one can define a probability density 
function with real variables. The master equation for this function is (to lowest order) 

BP BP 1 B 2 P 

"i-in-L^iy-^. (is) 

In this case, the Galilean boost becomes a symmetry of the space-time being considered. 
Going to a reference frame moving with the velocity V = R — L, one obtains 

dt 2^dxf l j 

The asymmetry due to driving has been removed. However, the asymmetry due to pushing 
remains intact since the boundary condition does not change under this transformation. 
The rest of this paper is devoted to the technical details and elaboration of the above results. 
In section 2 we prove that equations (0) and (|3]) actually give the above mentioned process. 
In section 3 we apply the coordinate Bethe ansatz and show that for this highly nonlocal 
process the iV-particle S matrix is still factorizable. We also find the integral representation 
of the iV-particle conditional probability distributions. In section 4 we show that in the two 
limiting cases /i = 1 and A = 1, the conditional probabilities can be expressed in closed form 
as N x N determinants. Section 5 is devoted to the mean field solution and discussion of 
the current density relation. In section 6 we calculate the drift and the diffusion rates for 
the two particle sector and compare our results with those of the ordinary partially ASEP 
[22]. Finally in section 7, we discuss the qualitative picture of the phases for open systems 
and also discuss the relation with ordinary ASEP in other updating schemes. 

2 The Master Equation and the Process 

The master equation is (||) supplemented by the boundary condition (|3|). For sectors of low 
number of particles (e.g. N=2,3), one can repeatedly use (|3|) to find the rates. For general 
sectors we use the following lemma. 
Lemma: The boundary condition @ implies 

P (x, x + 1, • • • , x + n — 1, x + n, x + n) = 

(1 — r n+ i)P(x, x + 1, ■ ■ • , x + n — l,x + n, x + n+ 1)+ 

r n+ iP(x — 1, x, • • • , x + n — 2, x + n — 1, x + n), (15) 

P (x, x, x + 1, ■ ■ • , x + n — 1, x + n) = 

(1 — l n+ i)P(x — 1, x, x + 1, ■ ■ • , x + n — 1, x + n, x + n)+ 

l n+ iP(x, x + l,x+ 2, ■ ■ ■ , x + n — 1, x + n, x + n + 1), (16) 
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where r n and l n are given in (|7|) and (||). 

The proof of this lemma is almost the same as that given in [23] and will not be repeated 
here. This lemma in fact states how to resolve the singularity in coordinates of a cluster of 
n adjacent particles with that of a single particle from right and left respectively. 
Consider now the master equation for a collection of n adjacent particles: 

d 



— P(x, x + 1, ■ ■ ■ , x + n — 1) 
at 



R 



n-l 



TiP{x — l,x, ■ • • ,x + i — 2, x + i — l,x + i + I, - • • ,x + n — 1) 

,i=0 
(n-l \ 

T{ J P{x, x + 1, • • • , x + n — 1) 



i=0 
n-l 



+ L 



£ L-i-iP(x, x + l,---,x + i — l,x + i + l,x + i + 2,---,x + n) 



u=o 



x, x + 1, • • • , x + n — 1) 



(17) 



It is now clear that the master equation (fj) and the boundary condition (H) imply the pro- 
cesses (ED-fBD. 



3 The Bethe Ansatz Solution 



3.1 The case of an infinite lattice 

The Bethe-ansatz solution to the master equation (§) is 

P(xt, ■■•,x N ;t) = e Et ^{xi, ■■■,x N ), 

where 

*(x 1 ,---,x N ) = '£A iT e ia <»>* 



(19) 



Here x and p denote n-tuples of coordinates and momenta, respectively, the summation runs 
over the elements of the permutation group, and A^s are coefficients to be determined from 



the boundary condition (H). Inserting (|T8|) in (0), we have 



N 

#£*(*!, 



N 



, Xj 



1, • • • , x N ) + L £ • • • , Xj + 1, • • • , x N ) = {N + (zi, • • • , x N ), 

(20) 



or 



(pj) 



.xn 



(21) 
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From this, one obtains (as one can remove a from the summations in the parenthesis) 

E = ^E( Pj ), (22) 

j 

where 

E( Pj ) = Re~ ip > + Le ip > - 1. (23) 

The next step is to determine the coefficients A a , so that the eigenfunction \1/ satisfies 
the boundary condition (£J). It is seen that nothing from the master equation enters this 
boundary condition. So the solution to this boundary condition is just what was found in 
[23], that is 

Ant = S{a(p i ),a(p i+1 ))A <J , (24) 

where a is an arbitrary element of the permutation group, and cr, is the generator which 
only interchanges pi and Pi+i- The elements of the two-particle scattering matrix, 

S jk := S( Pj ,p k ) = - i _ Ae . pj _ — k , (25) 

are thus sufficient to calculate the scattering matrix in the general iV-particle sector, and the 
latter factorizes in terms of the former. As the scattering matrix is just that obtained in [23], 
the same reasoning shows that there are no bound states. So we can write the conditional 
probability as 

P(x; t|y; 0) = / [ft ^e^-^*(x, p), (26) 

where the normalization of \l/ is chosen so that the coefficient of e ip x in \l/ is equal to unity. 
Integrations run over < pj < 2tt, and the poles of the scattering matrix are shifted from 
the integration region through the prescription 

S(pk,Pm) -»• S(p k + ie,p m ), k<m. (27) 
This ensures that in the physical region (xj < Xj + i , yj < Uj+x) we have 

P(x;0|y;0) = n<W (28) 
j 

Equation (26) is an integral representation for the conditional probabilities. More explicit 
information is obtained after calculating the integrals. For example in the two particle sector 
equation (|26| ) is written as 
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f d 2 p 

P(x 1 ,x 2 ;t\y 1 ,y 2 ;0) = J —e 



[E(p 1 )+E(p2)]t-i(piyi+p 2 y2) 



e i{vixi+p2X 2 ) t— e i(piX2+P2Xi) I (29) 

1 + e - Ae*w - fie-** 3 J ' 

Using the variables £ := e* Pl and r/ := e~ W2 , a contour integration yields 



P(x 1 ,x 2 ;t\yi,y 2 ;0) = e 2t J2 



n\(n + y 2 — x 2 )\m\{m + x\ — y x )\ 



\ p J n!(n + y 2 - x\ + q)\m\(m + x 2 - y\ +p)\ 
\(n + y 2 -Xi + q) //(m + x 2 - y x + p) 



1 



Lt Rt 



(30) 



Here all of the summations run from zero to infinity. A simple calculation shows that in the 
limit L = 0, the result of [23] is obtained. 

3.2 The case of a periodic lattice 

On an infinite lattice, the set of momenta of the eigenfunctions of the Hamiltonian are 
continuous. On a finite lattice, however, this set is discrete. To obtain this set, consider a 
lattice of M sites, on which N particles live. Now, another boundary condition should be 
added [25] 

V(xi, x 2 , ■ ■ ■ ,x N ) = ty(x 2 , x 3 , ■ ■ ■ , x N , Xi + M). (31) 

This means that one cannot unambiguously define the first particle: one can interpret the 
first particle as the last one, provided its coordinate is enhanced by M, the period of the 
lattice. Applying the boundary condition (]31) ) on the eigenfunction (0), we have 

Y^t4 e*[ <J (P 1 ) Zl +°"(P 2 ) Z2 ~ l ^{Pn)xn] _ e i[v(pi)'X2+'?(p 2 )x3-\ \-a(p N )(x!+M)] 

a a 

— ^2 e i[v(P2)x 2 +a(p-i)x 3 -{ \-cr(p 1 )(x 1 +M)} ^32) 



where 
This yields 
But, 



°o(Pi,P2,--- ,Pn) ■= (P2,P3, - ■ ■ ,Pn,Pi)- (33) 
A a = A uuo e lM ^\ (34) 
a = a x ---a N -x. (35) 
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So, using (Gl), 



4 4 

S1<J(TQ -^UUX ■-•CTJV— 1 

= A^o-cri • • • (Tiv_2(pjv-i), (rax--- (?n~2{pn)} ■ ■ ■ S[a(px), a(p 2 )] 

= A tJ S[a(px),a(p N )]S[a(px),a(p N „ 1 )]---S[a(px),a(p 2 )\. (36) 

Combining this with (|34[), we arrive at 

e -.AM Pl ) = 5[(T ( pi)) a{pN)] . . . S [a( Pl ), a(p 2 )\, (37) 

which can be written as 



e 



-iMp 



k 



US(p k , Pj ). (38) 

These are the Bethe equations of the system, the solution of which provides the allowed set 
of discrete momenta. Note that the driving parameter enters only the energy equation, as 
in (23), and the pushing parameter enters only the Bethe equations. Denoting e~ rpk by z^, 
eqs. ( p8|) can be rewritten as: 

4' = n r H + " z - " 1 < 39 > 

j+k 1 - AZ k - V Z j 

This system of equations have a symmetry, namely it is invariant under z — > z~ l , A — > fi. 
This means that if the set z a := {z%} are the quantized momenta for the (A,/i ) system, 
then the set u a := {{z^)~ 1 } are the quantized momenta for the (fi, A ) system. Thus if we 
know the spectrum of the former system (see eqs. (22) and (23)): 

E a = Y.(Rzf ] +L(zf ) y'-l), (40) 

3 

then the spectrum of the latter system is also known 

^ a =E(^j a) r 1 +^ a) )-i) (4i) 

i 

In particular this means that part of the analysis of the spectrum of ordinary ASEP (L = 
fi = 0) which has been done by Gwa and Sphon [25] , can be applied to the drop-push model 
[24]. 

4 Closed Form of the conditional probabilities in the 
limiting cases A = 1 and /i = 1 

In the special cases A = 0, 1, one can use a determinant ansatz for the conditional probabil- 
ities [22,23]: 

P(x;t|y;0) = e- 7Vi det[G(x;t|y;0)], (42) 



S 



where G is an N x N matrix with elements 



Gij(x] t\y; 0) = gi-j{xi - yf t). 



(43) 



Inserting (|^) in (g), one obtains 

d 

-Q t Gij{x; t) = R G^x - 1; t) + L G^x + 1; t). 

The equations obtained by the boundary condition (0) (for A = 0, 1) are 

( g k -i(x; t) = g k -i(x - 1; t) + (3g k (x; t), A = 
\ gk+i(x; t) = g k +i(x + l;t) + (3g k (x; t), A = 1. 

Writing (HI) in the form 



g k (x; t) = R g k (x - 1; t) + L g k (x + 1; t), 
and introducing the ^-transform 

g k (z,t) := J2 zX 9k(x;t), 

X 

we obtain from (f^J) the following 



(44) 



(45) 



(46) 



(47) 







-g k (z,t), A = 



9k-l(z,t) : _ 



(4S 



while (ES) yields 



^(z,t) = e ( &+L /^ fe (^,0). (49) 
We also note that g k (z,t) is the z-transform of the one-particle sector probability; so 



g o (z,0) = Y, zX Ko = 1. 
Combining; (|48[), (pE9"D, and flSOD, we arrive at 



(50) 



9h(z,t) = < 



D (Rz+L/z)t 



,(Rz+L/z)t | 



1 - z 



P 



A = 
A = 1. 



(51) 



V 1 - l / z ) 

The parameter /3 drops from the determinant, so that one can set it equal to an arbitrary 
number; we set it equal to unity. From this, one obtains 

k \ (Rt) m (Lt) n 



9k(x;t)= J2 



n—m-\-x | 



m,n — 



,n — m + x j m\n\ 



A = 0, 



(52) 
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and 

9k(x;t)= £ (-!)»—( "* V Lt rm n A = 1. (53) 
^ \n — m — x J mini 



m,n— 



Note that 

9 t^ L \ X ;t)= 9 it°' L ' R \-^t), (54) 

which is a special case of the symmetry under reflection. By this, we mean that the system 
of equations (2) and (3) is invariant under the following transformations: 

Xi > — Xjy+l-i, R ► L, A > [L 



5 Steady state of the system on a Ring 

In this section we consider a ring of N sites on which M particles are hopping. The steady 
state of this system is the one in which all the configurations have equal weights. Thus all 
the steady state probabilities P{x\, 22, ■■■Xm) are equal to a constant. Stationarity of this 
measure is proved by noting that P(xi, X2, ■■■X^) = constant, satisfies both the master 
equation (0) and the boundary condition (^). Uniqueness of the measure is ensured by 
connectivity of the process, i.e. the fact that every configuration can be reached from any 
other by a sequence of transitions [26]. In this state one can calculate all the correlation 
functions by simple combinatorics. If n k is the random variable at site k which is 1 if it is 
occupied and if it is vacant, then it is well known [15] that 

M M(M-l) M(M-l)(M-2) 

<"i>=^, >= N(N _ iy <n j n k n l >= N{N _ 1){N _ 2 y etc. 

(55) 

In the thermodynamic limit, when M —>■ 00 and iV — > 00 with M/N = p, this steady state 
approaches the uncorrelated steady state given by the mean field solution. What we want 
to do in this section is to find the current density relation for such a steady state. In general 
one can find the equation for the rate of change of the average density at site k, either by 
going to the Hamiltonian formalism and using the equations -g < n k >=< [n k , H] >, or by 
just looking at the process and determining the various ways in which this density decreases 
or increases. We follow this second approach which is more transparent and intuitive. The 
rate of change of density of particles can also be written as a continuity equation, i.e. ^ < 
n^ >= Jk-i — Jki where J k is the current through site k. This current is the algebraic sum 
of a positive and a negative current 

J k := RJ£ - LJj^, (56) 

where J~£ and J^T are due to the hopping of particles to the right and to the left respectively. 
Due to the pushing effect, both of these currents are non-local. The explicit expression of 



10 



Jjt and J k are: 



j£ = r < n k (l - n k+ x) > +n f < n fc _in fe (l - n fe+ i) > + < n fc n fe+ i(l - n fc+2 ) > 
+ r 2 f < n k _ 2 n k -in k {l - n k+1 ) > + < n fe _ 1 n fc n fe+1 (l - n k+2 ) > 
+ < n fe n fe+ in fe+2 (l - n k+3 ) > J H , (57) 



and 



J; 



= l < (1 - n fc )n fc+ i > +ij < (1 - n fe )n fc+1 n fe+2 > + < (1 - n k -i)n k n k+x > ^ 
+ fef < (1 - n k )n k+1 n k+2 n k+3 > + < (1 - n k ^i)n k n k n k+ in k+2 > 
+ < (1 - n k _ 2 n k _ x n k n k + x ) > J H . (58) 

A typical term like < n k _in k (l —n k+ i) > in the expression of JjT, is in fact the probability of 
the configuration 110 on sites k — 1, k, and k + 1, respectively. We know from (^) and (0) that 
this configuration changes with rate ri to Oil on the same sites, hence this term contributes 
to the current Jjt. A similar interpretation is true for other term < n k n k+ i(l —n k+2 ) >, etc. 
In the uncorrelated steady state, the above currents are calculated to be 



J + :=P(1"P)E 



(n + l)p r 



n=0l + H 



and 



from which one obtains: 



and 



J- 



J := P(l " P) E 



(n + l)p r 



+(*)' 



(59) 



(60) 



P(l-p)(l + yP) + o(f) 2 , 



P- 



A 

— / 

P 



(61) 



^a + -p( P -2)) + o(A)2, 



T^(l + ^(p-2)) + o(f) 2 , 



P 



P 



2A 



»1, 



« 1 



£ = 1, 



(62) 



A\2 



p(l-p)(l + — p) + o(±) 



Consider J + : It is seen that as far as the pushing effect to the right is small (f << 1), the 
standard mean field current of the ASEP gets only corrections of the order f . For medium 
pushing, when f ~ 1, the current is exactly equal to the density and the interesting point is 



11 



that contrary to the case of ASEP, even when the lattice is filled with particles (p = 1), there 
is a nonzero current due to pushing. At very strong pushing » 1), the current even 
diverges when the lattice is filled . It is now instructive to consider only the leading terms 
of the total current Jk in different regimes. In the steady state, this current is independent 
of the site number k and hence is denoted by J. We find from ( |B~TD and (|62|) the following 



Remark: Note that our model incorporates only the totally ASEP as a special case . There- 
fore when L ^ 0, the above results should not coincide with those of the partially ASEP in 
the limiting case \x = 0. 

Consider the case j << 1, where we have strong pushing to the left and weak pushing to 
the right. It is seen that for L > R the current is always negative, which is expectable on 
physical grounds. However, for L < R, when the driving force is to the right, the two effects 
act in opposite directions. The current is positive as long as p < p c := 1 — J~^. At p = p c 
the current vanishes and for p > p c , the pushing effect takes over and the current becomes 
negative. 



6 Drift and Diffusion Rates in the Two Particle Sectors 

In this section we want to study the interplay of driving and pushing in the behaviour of two 
important dynamical quantities, namely the drift and the diffusion rates. More generally, 
we study the long time behaviour of the quantities 4 < X > and 4(< X 2 > — < X > 2 ). 
Our starting point is the exact calculation of conditional probability of the the two particles 
being a distance x apart, given their initial separation y. Denoted by P r (x; t\y; 0), it is given 




(63) 



as: 



oo 



P r (x;t\y,;0) 



P(x 2 - x,x 2 ;t\0,y;0) 




+ S(p!,p 2 )e 



-IP2X 




(64) 



where we have used: 



—ip 



(65) 
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From these we arrive at 



P r {x; t\y; 0) = e~ 2t [l y - x (2t) + l y+x -i(2t)} , (66) 
where I n is the modified Bessel function of order n with the integral representation 

= I 7^e inp+ucosp . (67) 

It is interesting to note that this probability is independent of the asymmetry- and drift- 
parameters. Another way to see this result, is to derive the equation of evolution for P r . 
To do so, begin from the master equation (§) for two particles. Using the definition of 
P r {x; t\y; 0) now abbreviated to P r {x), we arrive at 

P r (x) = P r (x - 1) + P r (x + 1) - 2P r (x). (68) 

The boundary condition (^) is transformed into 

P r (0) = P r (l). (69) 

It is seen that the driving and the pushing parameters are absent in this equation. The 
physical explanation for the absence of /x or A is that, when two particles push each other 
they do not change their inter-particle distance. This distance increases by one unit with 
rate R + L (particle 2 hoping to the right or particle 1 to the left), and decreases by one 
unit with the same rate (particle 2 hoping to the left or particle 1 to the right), and since 
R+ L has been rescaled to unity, the driving parameters do not appear in these equations 
either. We should stress that this is not the case in more than two-particle sectors and the 
probabilities for relative distance in these sectors do indeed depend on the above parameters. 
The next quantities we calculate are the average velocities of particle 1 and particle 2. Note 
that particles can not overtake each other and that they keep their initial order at all times. 
We have: 

■ •: Y^xl'i.r) (70) 

X 

where P{(x) is the probability of finding particle i at site x . The master equation for these 
probabilities are obtained from (0) and (0), using the definitions: 

oo 

P 1 (x):= ]T P(x,x 2 ), (71) 

x 2 =x+l 

and 

P 2 (x):= J2 P(xi,x). (72) 

x±=— oo 

This calculation finally leads to the following equations 

Pi(x) = R{P 1 (x-l)-P(x-l,x)+fxP{x-l,x)}-(x — >x + l) 
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+L{Pi(£ + 1) + \P(x + 1, x + 2)} - (x — > x - 1), (73) 

and 

P 2 ( x ) = P{P 2 (x - 1) +/iP(x - 2,x - 1)} - (x — >x + l) 

+L{P 2 (x + 1) - P(x, x + 1) + AP(x,x + 1)} - {x — >x- 1). (74) 

We have written these equations in this unsimplified form in order to convey their simple 
physical meaning. In fact they can also be obtained by intuitive reasoning . Consider for 
example eq. (|73|). The first two terms in the curly bracket of the first line are due to particle 
1 at site x — 1 hopping to an already vacant site at x and the third term is due to particle 
1 at site x — 1 hopping to site x and pushing the already present particle 2 at this site to 



the right. Other terms have similar meaning. One now obtains from (70) and (^,0) the 
following 

^^ = P-L-AP r (l), (75) 

and 

d< * t 2> =R-L + »P r (l). (76) 

One can even derive these equations from the beginning by physical reasoning without using 
the master equation. For example we know that the hopping rate of particle 1 to the right is 
normally R unless it is one site behind particle 2 where its hopping rate will be Rfi. These 
two terms can be combined to give a positive contribution R + (P/i — P)P r (l) to the average 
velocity of particle 1. On the other hand its hopping rate to the left is normally L unless 
particle 2 is exactly one site to its right, where its hopping rate becomes L + LA. These 
two terms are then combined to give a negative contribution —L — LXP r (l)) to the average 



velocity. Adding these, one obtains (75). The same kind of reasoning can give (76). Using 
the definitions 

< r >:=< x 2 > - < xi > < X >:= i(< x 1 > + < x 2 >) (77) 



we find 



and 



d < r > 



= H-L + ^^>. (79) 



From these, one obtains 



< X >=< X > +(R - L)t + ^— — (< r > - < r > ), (80) 
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where the subscript refers to initial conditions. From (|66|), and using the asymptotic 
behaviour of the modified Bessel functions, one can obtain the asymptotic behaviour of 
these expectation values. We have 

1 



P r .(l) = ^= + 0(r 3 / 2 ). (81) 



Then, we obtain 

<r >=C + 2^ + 0(r 1/2 ), (82) 
where C is a constant depending on the initial conditions. So, 

< X >=< X > +(R - L)t + ^— ^ ( C + 24/-- < r > ] + 0(r^ 2 ). (83) 
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= R-L + ^ + 0(r 3 / 2 ). (84) 



At long times we have: 



It is seen that to leading order the drift rate defined as V := 4 < X > is only controlled by 
driving, and pushing has only a sub-leading effect. To calculate the diffusion rate we proceed 
as follows: 

where 

< Xi ? >= ^a; 2 ^(a;), i = l,2. (86) 

a; 

Using ( [73] ) and ([71]), we arrive at 

d < x 2 > 
dt 

and 



l + 2(i?-L) <xi >+A(L-i2)P r (l)-2A5^a:P(x,z + l), (87) 



d < a; 2 . > 



1 + 2(R - L) < x 2 > +fi(R - L)P r (l) + 2fiY,x P(x, x+1). (88) 



dt 

From these we obtain 

^^'V^^ = l-(,-A)<X>P r (l) + 1 + ("-f fi - L ) p r ( 1 ) 

+ ^^E(2i-1)-P(i-1,i) (89) 

The value of the quantity J2 X (^ X ~ l)P(x — 1, x) is calculated in the appendix . Inserting its 
value in the above formula, using the asymptotic form of the modified Bessel function, one 
obtains 

1: _ _d(<X*>-<X>*) , „ . W l 1 
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We see that pushing has a leading effect on the diffusion rate. The physical explanation 
behind this is that although the average distance between the particles grows with time 
as t~2 (see (p2|)), the width of the probability distribution (wave packet) of each particle 
increases also with rate £2, thus the wave packets always overlap and there is always a finite 
probability that the particles push each other. 

7 Discussion and Summary 

We have defined a generalized asymmetric exclusion process with random sequential up- 
dating in which particles besides hoping randomly to the left and right can also push their 
neighboring particles with rates depending on the number of these particles. We have shown 
that this model, although governed by a very non-local Hamiltonian of the spin chain type, 
is exactly solvable via the coordinate Bethe ansatz. (For the type of the Hamiltonian see 
[23]). Throughout the paper we have tried to study the interplay between the two sources of 
asymmetry, the one due to driving and the one due to pushing. Due to the pushing effect, in 
our model, cluster of particles can also hop to right and left. This is similar to what happens 
in sequential updating schemes [20,21], where in one complete updating of the lattice, clus- 
ters of particles move. However, we remark that here we have this effect in continuous time 
and not discrete time. Specially, in sequential updating scheme, when the total updating 
operator is the product of local updating operators, the probability of hopping of clusters in 
one complete update, turns out to be a power of the hopping rates of single particles while 
in our model this is not so. That is why our Hamiltonian is very non-local. 
The fact that despite this non-locality the system has an exact solution and its S-matrix is 
factorizable is interesting. What remains to be done for this model is to study its steady 
state ( particularly in the totally asymmetric case L=0) and its phase structure on open 
systems when particles are injected and extracted at the open ends, to see how the simple 
phase diagram of the ASEP will be modified due to pushing. Again, due to the non-locality 
of the Hamiltonian of the process, the conventional technique of Matrix Product Ansatz can 
not be applied to this problem. Qualitative pictures may be obtained along the work of ref. 
[27]. However, one should first decide as to how to add boundary terms to this process, that 
is if particle arrive only at the already vacant boundary site or else, they can also push a 
cluster of particles of arbitrary size, already present there. If this is so, then there will be 
no longer a genuine difference between the boundary terms and the bulk. We believe that, 
due to these complications and the nonlocal character of the process, this problem deserves 
a separate study. 
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8 Appendix 

To calculate the quantity J2x(2x ~ l)P(x — l,x) needed in section 6, we proceed as follows: 



^2(2x-l)P(x-l,x) = J 



d 2 p 
An 2 i 



,Et-ip-y 



d d 
+ 



E 



D i(pi+P2)x-ipi 



+Su \dpi + dp 2 ) ? 6 

/id 2 p 



dpi dp 2/ x 

i(pi+P2)x-ip2 



e -ipi ( _|_ _9_) e Et-ip-y 

\dpi dp 2 



{dpi dp 2 . 



(S 12e Et-i P .y 



e- 2t {( yi + y 2 )[l y (2t) +l y -i(2t)} 

+t(R - L)[l y+1 (2t) + I„(2t) + l y -i(2t) + l y - 2 (2t)} 

oo 

(20]}. (91) 



fc=0 



Using the identity 



we arrive at 



E Wt) = e 2t , 



(92) 



J2C2x-l)P(x-l,x) = e- 2t {(y 1 + y 2 )[l y (2t) + I y _i(2t)] 

+t(i? - L)[I, +1 (2t) + 1,(20 + Vi(20 + l y -2(2t)}} 
v-i 

+ (/i - A){l - e- 2t £ [In(20 + W(2t)]}. (93) 

n=0 
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